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Abstract 

Weekly coordinate time series for the all EPN (European Permanent GPS Network) stations 
for the period of observations beginning from GPS week 834 was obtained by a reprocessing 
of the official EPN solutions. Comparison of this series with EPN Central Bureau (EPN CB) 
and the previous lAA two-year solution for selected EPN stations obtained by reprocessing of 
the original observations is performed. Results of comparison show that new solution based on 
fiducial-free strategy is most likely free of seasonal errors. Now the new lAA EPN solution is 
being computed regularly (on availability of the official EPN solution) and is available to any 
interested group. 

1 Introduction 

GPS observations collected from the European GPS network are widely used for geodesy and 
geodynamic researches in the region. The network is coordinated by the EPN CB which also 
provides the official analysis of the observations. Two EPN weekly solutions are available. The 
first one is computed at the Bundesamt fiir Kartographie und Geodasie, Germany (previously 
at the Center for Orbit Determination in Europe, Astronomical Institute of the University of 
Bern, Switzerland) and distributed as SINEX files. Hereafter this solution is referred to as EUR. 
Unfortunately, this solution is not suitable for high-accuracy geodesy applications because the using 
of fiducial approach and periodic change of the reference system cause jumps in coordinates and 
distortion of the network. Direct use of this solution can lead to some confusions (see e.g. ^). 
Besides the EUR solution does not contain any information about displacement of the fiducial 
stations. 

The second solution is computed at the EPN CB by reprocessing the previous one and seems 
to provide high-quality information about the movement of the all EPN stations. Unfortunately, 
this solution is not distributed in SINEX files, and it is difficult to use it in scientific analysis. 

For these reasons, several years ago the lAA undertook a special project aimed at the compu- 
tation of an independent coordinate time series for all the EPN stations. A detailed description of 
the project was given in [2]. 

At the first stage the original GPS observations for a selected EPN subnetwork at two-year 
interval were reprocessed using fiducial-free strategy. Basic theoretical background for this approach 
can be found in [HllSIEliniEilHlin]- Previously this strategy was tested during processing of two 
Baltic Sea Level campaigns E31 Obtained coordinate time series appeared to be more stable in 
sense of random and systematic errors [SUH]. It is hereafter referred to as II. 

Unfortunately, this way of reprocessing requires too much resources and we tried two another 
approaches. The first one is based on reprocessing of existing EUR solutions (the second stage of 
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the lAA project). The processing strategy is described below. Obtained coordinate time series was 
compared with other EPN solutions and the previous lAA solution II. Results of comparison show 
that obtained solution is most likely free of seasonal errors. 

At the third stage of the our project we plan to obtain an EPN coordinate time series by an 
independent combination of individual solutions provided by the EPN Analysis Centers, also using 
fiducial-free approach. This work is under development and is planned to be completed by the end 
of 2003. 

2 Processing strategy 

Our processing was made in two steps. At the first one the EUR solutions are de-constrained using 
the a priori coordinates and covariance matrices contained in the SINEX files following the strategy 
proposed in [TT] . 

After this the transformation of the obtained free network solution to the ITRF2000 is made. 
Unlike EPN CB strategy, we use for the transformation all the stations presented in the solution 
using weights dependent on their position accuracy (taking into account both errors in position 
and velocity) in the EUR solution and in the reference ITRF catalogue. 

Two solutions were computed: with 6 and 7-parameter Helmert transformation. Solution ob- 
tained with 6-parameter transformation is hereafter referred to as 12, solution obtained with 7- 
parameter transformation is referred to as 13. 

Time series of transformation parameters are shown in Figure^ One can clearly see the seasonal 
and other peculiarities in these time series which should be investigated separately. 

Analysis of the solutions obtained by this method showed that the obtained coordinate time 
series have a similar quality with the EPN CB solution with some small discrepancies which can be 
explained by details of used approaches. However analysis of the errors in the station coordinates 
reveals very large irregularities, caused most probably by inconsistency in SINEX blocks (e.g. wrong 
scaling of covariance matrices). 

To reduce the inconsistency of coordinate errors we used a re-scaling of computed covariance 
matrices. The re-scaling factor was chosen in such a way that the mean coordinate error of non- 
fiducial stations in the EUR solution is equal to the mean coordinate error of the same stations of 
the new solution. The plots of mean coordinate errors of our solution with and without re-scaling 
of covariance matrices are shown in Figure |21 

3 Comparison and conclusions 

Six-year coordinate time series for all EPN stations were computed and compared with two other 
solutions (EUR solution and our previous two-year solution for selected european stations). Several 
examples are presented in Figure |21 It is seen that all series provide determination of main details 
in behavior of station position. 

A comparison of four solutions mentioned above is presented in Table ^ The table contains 
results of determination of week-to-week repeatability (Allan variance) interpreted as random error 
and amplitude of seasonal term in variation of station coordinates in the local ENU system. 

Comparison shows that 13 solution obtained using 7-parameter Helmert transformation to the 
ITRF2000 provides minimum random error and seasonal variations. Of course, the latter may mean 
merely loss of geophysical signal, but comparison with EPN CB series (available at the EPN Web 
site) and global solutions (T. Springer, private communication) shows that most likely seasonal 
terms observed in the EPN solution is caused by systematic errors induced by errors in modeling 
of position of fiducial stations. 
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Figure 1: Helmert parameters time series for transformation w.r.t. ITRF2000: with 6 parameters 
(left) and with 7 parameters (right) (R — rotation parameters, T — translation parameters, D — 
scale); vertical bars on the plots present parameters errors. 
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Table 1: Statistics for 10 stations presented in all solutions. 
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Figure 2: Mean coordinate errors without (top) and with (bottom) re-scahng, mm; 

Evidently, more thorough consideration should be made to make a choice between using 6 or 7- 
paramctcr Hclmcrt transformation of free network solution to ITRF. The first impression is that it 
is more reasonable to apply 6-parameter transformation to a global network, whereas 7-parameter 
transformation is more adequate to regional data. 

Now the lAA EPN solutions are being computed regularly using described strategy, on avail- 
ability of the official EPN solutions, and the results are available to any interested group (on request 
at the moment, later it will be put to the Internet). 

4 Aknowledgments 

Authors are very grateful to Matthias Becker and Daniel Ineichen for their help in understanding 
of strategy of computation of the EPN combined solutions. 



5 



EUREF 

WTZR, Up component 



40 
20 



-20 
-40 



850 900 950 1 OOO 1 050 1 1 OO 



WTZR_14201H010 



cfP" c,TS^ c^-^ q'f' q'S' ^fP" q'P q"-^ q«R ^-^l- ,TS- 

IAA-1 

WTZR, Up component 



40 
20 



-20 
-40 



40 

20 



-20 
-40 



850 900 950 1 0OO 1 050 1 1 0O 

IAA-2 

WTZR, Up component 



850 900 950 1 OOO 1 050 1 1 OO 



EUREF 

HERS, East component correcteci for NNR-NUVELI A 



40 
20 
O 

-20 
-40 



S50 900 950 1000 1 050 1 1 0O 



HERSj32l2H0a7 







: * Ea.'?l-a)m|»n&nl 





^'i" 



■.^ -.^^ vl" # 



j# ,0* a'^' .v^ 



IAA-1 

ERS, East component corrected for NNR NUVELI A 



40 
20 



20 

-40 



40 
20 
O 
-20 

-40 



850 900 950 1 OOO 1 050 11 OO 

IAA-2 

= RS, East component corrected for NNR-NUVEL1 A 



850 900 950 1000 1 050 1 1 OO 







EUREF 

ZWEN, Up component 


40 

20 

-20 
-40 


■.Y'hvr^ n,,r^ 






S50 


900 950 1000 1 050 1100 

ZW&N 12330HDD1 

.ffi- .od^ ^ 




ID 


















SI ^^P q<fl qTSi qfP- ^ ^jlP ^j'fP ^flP -.iSJ- v 

IAA-1 

ZWEN, Up component 




40 
20 


-20 
-40 








850 


900 950 1 000 1 050 1 1 00 

IAA-2 

ZWEN, Up component 


40 

20 


-20 
-40 








S50 


900 950 1000 1 050 1100 








EUREF 

SVTL, Up component 


40 

20 


-20 
-40 








850 


900 950 1000 1 050 1100 

SVTL 12330 MD01 

.fii- .ff^ .ofi^ ^ 






















^ '=1^ ^f*^ ^ tf^ ^ 

IAA-1 

SVTL, Up component 




40 

20 


-20 

-40 








S50 


900 950 1000 1 050 1100 

IAA-2 

SVTL, Up component 


40 

20 


-20 
-40 








S50 


900 950 1000 1 050 1100 



Figure 3: Some examples of the coordinate time series, mm (top to down: EUR solution, EUREF 
CB solution (copied from the web page tl2|), II solution, and 12 solution.). 
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